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III.  Method  of  Analysis 


Computer  Program 

In  order  to  properly  account  for  the  complex  response 
of  a  metal  subjected  to  high  temperature  cyclic  loading,  a 
specialized  finite  element  program  has  been  used.  VISCO, 
utilizing  the  Bodner-Partom  constitutive  law,  is  a  two- 
dimensional  viscoplastic  finite  element  program  capable  of 
accurately  modeling  cyclic  loading,  cyclic  hardening  or 
softening,  cyclic  relaxation,  and  cyclic  creep  [13]  .  It 
uses  constant  strain  triangles  for  both  plane  strain  and 
plane  stress  solutions.  VISCO  was  originally  developed  by 
Hinnerichs  [6,  7,  8],  with  accuracy  verified  by  Smail  [9], 
Keck  [10] ,  and  Wilson  [11]  . 

The  Bodner-Partom  viscoplastic  constitutive  equations 
in  VISCO  are  economically  solved  using  the  Gauss-Seidel 
indirect  solution  procedure  incorporating  an  optimum 
overrelaxation  factor  to  speed  convergence  [14].  An  Euler 
extrapolation  scheme  is  employed  for  the  numberical  time 
integration  of  the  Bodner  equations  [6] .  The  incremental 
solution  process  is  based  on  incrementing  time  directly. 
Loads,  strains,  stresses,  etc.  are  thus  incremented 
indirectly  and  are  solved  for  using  the  residual  force 
method.  The  time-step  size  is  checked  by  comparing  results 
with  prescribed  constraints.  These  constraints  are  the 
allowable  amounts  of  change  in  stress,  o  ,  and  strain, 

1  U1j 
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Table  2 . 1 


Bodner  Coefficients  for  IN-100  at  1350°  F 


Material 

Parameter 

Description 

Value 

E 

Elastic  modulus 

26 . 3x102KSI 
( 18 . 133xlO^MPa) 

n 

Strain  rate  exponent 

0.7 

D 

o 

Limiting  value  of 
strain  rate 

,  „4 

10  sec. 

Z 

o 

Limiting  value  of 
hardness 

915 .0KSI 
(6304  MPa) 

Z1 

Maximum  value  of 
hardness 

1015. 0KSI 
(6993  MPa) 

Z2 

Minimum  value  of 
hardness 

600. 0KSI 
(4134  MPa) 

m 

Hardening  rate 
exponent 

2.57KSI-1 
(.37273  MPa  ) 

A 

Hardening  recovery 
coefficient 

1 .9xJ0~3 
sec 

r 

Hardening  recovery 
exponent 

2.66 

(1  KBAR  =  100  MPa  =  14 

.  5 0 4KS I ) 

A  and  r 


is  the  stable  non-work  hardened  value  of  Z. 
are  material  constants  for  IN-100  at  1350°  F  [12].  All 
the  Bodner-Partom  material  constants  are  listed  in  Table 
2.1.  Note,  these  values  are  constant  only  for  the 
specified  temperature  and  material. 

It  is  pointed  out  that,  for  fatigue  type  loading,  the 
load  can  change  sign  and  effect  the  normal  stress  strain 
relationship.  The  assumption  made  here  is  that  the 
material  is  isotropic.  Consequently,  the  properties  pre¬ 
sented  for  a  tension  stress  strain  function  are  duplicated 
when  the  specimen  is  stressed  into  the  compression  range. 


strain  rate,  Z  is  the  measure  of  material  hardening,  and 
exponent  n  is  a  rate  sensitivity  parameter.  Only  Z  depends 
on  the  deformation  history  of  the  material.  Z  is  assumed 
to  be  a  function  of  plastic  work,  W  ,  such  that 


mW 


Z  =  Z1  +  (Z  -Zj)  EXP(- 

o 


:  2 . 1 9 ) 


Z ^  is  the  maximum  expected  value  of  Z,  ZQ  is  the 
initial  value  of  Z  correspondning  to  the  reference  point 
from  which  plastic  work  is  measured,  and  m  is  a  material 
constant  that  controls  the  rate  of  work  hardening.  In  this 


thesis,  W  the  plastic  work  done  relative  to  some  initial 
P 


state,  is  the  only  variable  in  Equation  (2.19)  [6],  At  low 

temperature 


W  =  IS.  .  eijP(lt  (2.20) 

P  ij  J 

However,  this  analysis,  based  on  high  temperature 
conditions,  requires  a  thermal  recovery  term  which 
effectively  reduces  the  increase  of  material  hardening  due 
to  plastic  deformation: 


W 

P 


e .  .^dt  + 
ij 


Z  rec  dt 
mCZj-Z) 


(2.21) 


where 


I 


Bodner-Partom  Constitutive  Law 

The  Bodner-Partom  constitutive  law  is  based  on 
dislocation  dynamics  which  suggests  a  continuous  flow 
relationship  exists  between  stress  and  viscoplastic  strain 
starting  at  the  onset  of  loading  [5].  As  such,  the 
Bodner-Partom  equations  are  able  to  represent  the  principal 
features  of  cyclic  loading  behavior,  including  recovery 
upon  stress  reversal,  cyclic  hardening  or  softening,  cyclic 
relaxation  and  cyclic  creep  [6,  13].  The  formulation  of 
the  equations  is  arrived  at  by  squaring  the  Prandtl-Reuss 
relation  (2.13) 


e.  .p  e.  ?  =  \  s. .  s. . 

elJ  elJ  A  1J  1J 


[2.14) 


now  substitute  into  this  equation  the  following 


J2  *  t  Sij  SiJ 


,p  -i  e.  P£.  ,P 

}-  2  1J  1J 


yields 


D2P  =  ^  J2 


(2.15) 

(2.16) 

(2.17) 


whe 


P  . 


re  D2  is  defined  as  the  second  invariant  of  plastic 


strain  rate.  Bodner  and  Partom  expressed  D2F  as 

2  nvn  r  /•?  N  /  n+1  \  1  (2.18) 


D2p  -  D/  EXP  (f)l 


This  expression  is  based  on  extensive  experimental  data  and 
has  been  modified  to  fit  results  found  by  several 
researchers  [12].  Dq  is  the  limiting  value  of  plastic 
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u 


d 


—  j 

2G  J2 


(2.9) 


G  is  the  shear  modulus.  Under  the  Von  Mises  distortion 
energy  theory,  yielding  begins  when  the  distortion  energy 
in  the  multiaxial  case  equals  the  distortion  energy  at 
yielding  in  the  uniaxial  case.  In  terms  of  principle 
stress , 

J2  =  I  t (cf1-o2)2  +  Ca2-°3)2  +  (a3-°l)2}  (2.10) 


at  the  yield  point  in  uniaxial  tension,  J2  reduces  to 

J,  4o  2  (2.11) 

2  3  ys 

where  o yg  is  the  uniaxial  yield  stress.  Thus,  yielding 
occurs  when 


\  [(oro2)2  +  (o2-o3)2  +  (a3-o1)2]  =  ays2  (2.12) 

Now  we  have  a  criterion  to  predict  when  yielding 
begins;  however,  we  still  need  a  way  to  relate  viscoplastic 
strain  to  stress  in  a  manner  that  conveys  the  flow  nature 
of  the  problem.  This  requirement  is  satisfied  by  the 
Prandtl-Reuss  equations  in  rate  form  [15]. 

e.  ,P  =  AS..  (2.13) 

ij  iJ 

where  e.  ,p  are  the  components  of  the  deviatoric  viscoplastic 
ij 

strain  tensor,  and  are  the  components  of  the  devia¬ 

toric  stress  tensor  and  A  is  a  positive  proportionality 
constant.  Eqn.  (2.13)  will  be  used  in  subsequent  develop¬ 
ment  . 
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physical  properties  of  a  coldworked  metal.  Recovery 
reduces  the  effect  of  strain  hardening,  the  process  whereby 
the  stress  required  to  produce  further  plastic  deformation 
is  increased  because  of  prior  plastic  strain.  Creep  is  the 
time  dependent  strain  that  occurs  even  when  stresses  well 
below  the  yield  stress  exists.  If  a  specimen  is  deformed 
and  then  held  i;  a  fixed  position,  the  stresses  will 
gradually  decrease  due  to  the  creep  process;  this  is  called 
relaxation  [5] . 

A  model  that  includes  purely  elastic  as  well  as 
viscoplastic  behavior  is  called  elasto-viscoplastic  [16]. 
Total  strain  in  this  model  separates  elastic  (reversable) 
strain  from  viscoplastic  (irreversible)  strain: 


e  .  . 
ij 


(2.1) 


By  taking  the  time  derivative  of  equation  (2.1),  an 
expression  for  the  total  strain  rate  is  obtained. 


e.  . 

ij 


(2.2) 


The  elastic  strain  rate,  is  related  to  the  stress 

rate  simply  through  the  time  derivative  of  Hookes  law. 
Unfortunately,  viscoplastic  strain  rate  must  be  related  to 
stress  by  some  other  means. 

Conventional  plastic  deformation  begins  at  yield,  and 


is  dependent  on  the  yield  criterion.  Many  theories  are 
available  to  define  the  stress  level  at  which  plastic 


I I .  Viscoplastic  Theory 


Viscoplasticity  combines  two  inelastic  strain  groups, 
conventional  plasticity  and  time-dependent  deformation. 

The  conventional  theory  of  plasticity  is  characterized  by 
an  irreversible  strain  which  is  not  time  dependent  and 
which  can  only  be  sustained  once  yield  stress  has  been 
reached  [16] .  This  theory  has  been  satisfactory  since  time 
rate  effects  are  generally  not  important.  However,  metals, 
especially  at  high  temperatures,  exhibit  significant  time 
dependent  deformation  which,  like  conventional  plasticity 
is  irreversible.  These  two  permanent  strains  must  be 
combined  in  a  unified  viscoplastic  model  to  adequately 
study  high  temperature,  time  dependent  problems. 

A  discussion  of  viscoplasticity  must  include  the 
effects  of  load  rate  and  temperature.  In  most  metals,  high 
rates  of  load  application  result  in  less  plastic  flow  and 
higher  stress  fields  than  low  load  rates.  This  is  called 
rate  sensitivity  [17].  The  magnitude  of  this  effect  is 
relatively  minor  at  room  temperature  but  increases  rapidly 
as  temperature  increases.  This  suggests  the  analogy  of 
viscous  flow  in  a  fluid. 

Three  time-dependent  processes  also  increase  in 
importance  as  temperature  rises:  recovery,  creep,  and 
relaxation.  Recovery  is  defined  as  the  restoration  of  the 
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experimentally  determined  by  Stouffer  [12].  The 
Bodner-Partom  viscoplastic  flow  law  has  the  capability  to 
predict  the  behavior  produced  by  cyclic  effects  [13]  .  This 
flow  law  is  integrated  through  time  by  an  Euler 
extrapolation  scheme  [14] ,  and  the  law  is  incorporated  into 
the  finite  element  program  by  utilizing  the  residual  force 
technique  [15] . 

Finite  element  models  of  a  compact  tension  specimen 
used  by  Wilson  [11]  and  a  center  crack  specimen  used  by 
Hinnerichs  [6]  were  investigated  in  order  to  compare  the 
effect  of  different  load  geometries.  Load  was  input  as  a 
sawtoothed  stress-time  pattern  of  constant  amplitude  with 
zero  mean  load  (load  ratio  or  "R-ratio”  of  -1.0).  For 
comparison,  the  compact  tension  specimen  was  studied  at  two 
cyclic  load  frequencies,  2.5  Hz.  and  .167  Hz.  Maximum  load 
amplitudes  for  each  case  was  set  to  provide  a  stress 
intensity  factor  (Kj)  at  the  crack  tip  of  35  ksi/in.  The 
study  was  then  repeated  at  a  Kj  of  45  ksi/ in.  Hopefully, 
this  analysis  will  contribute  to  the  understanding  of  the 
failure  mechanism  in  critical  high  temperature  jet  engine 
components  and  help  make  a  retirement-f or-cause  program 


more  effective 


cycle  predictions  could  be  made  for  components  with 
subcritical  flaws.  Only  those  components  with  a 
quantifiable  critical  flaw  size  would  therefore,  need 
to  be  retired  [ 5 ) . 

Approach 

Due  to  requirements  imposed  by  a  retirement-for-cause 
program  for  engine  components,  a  significant  volume  of 
research  has  been  devoted  to  high  temperature  fracture 
mechanics  and  the  technology  appears  to  be  maturing 
rapidly  [1] .  The  study  of  material  behavior  has  evolved 
from  high  temperature  creep  analysis  [6,  7,  8,  9]  to  high 
temperature  elasto-VISCO  plastic  analysis  under  cyclic 
loading  [10,  11]. 

This  new  study  compares  the  effects  on  material 
behavior  resulting  from  different  cyclic  load  geometries, 
frequencies,  and  amplitudes.  The  subject  material  is 
IN-100,  a  superplastical ly  forged  nickel-based  superalloy 
used  in  turbine  disks  for  the  F-100  jet  engine.  VISCO,  an 
in-house  computer  program  developed  by  Hinnerichs  [6]  was 
used  for  this  analysis.  VISCO  is  a  finite  element  program 
which  uses  constant  strain  triangular  elements  and  has  the 
capability  to  run  cyclic  loads. 

The  Bodner-Partom  flow  law  subroutine  in  VISCO  was  used 


to  model  the  plastic  flow  during  the  load  cycling.  The 
Bodner  material  parameters  for  IN-100  at  1350  F  were 


turbine  disk  containing  subcritical  flaws,  many  parts  that 
are  currently  retired  could  be  kept  in  service.  Parts 
would  be  inspected  at  intervals  determined  by  the  ability 
to  reliably  detect  flaws  and  provide  assurance  of  adequate 
safe  life,  and  thus  parts  would  be  returned  to  service 
until  they  could  no  longer  pass  inspection  requirements. 

The  Air  Force  calls  such  a  program  retirement-f or-cause  and 
is  currently  studying  its  feasibility  for  many  high  cost 
components  limited  by  low  cycle  fatigue,  like  F-100  jet 
engine  turbine  disks.  If  we  assume  a  15-year  engine 
lifetime,  a  retirement  for  cause  program  could  result  in 
F-100  engine  life  cost  savings  of  about  $249  million  [2]. 
Actual  testing  of  retired  turbine  disks  has  verified  that 
significant  cost  savings  can  be  safely  realized  using 
retirement-f or-cause  procedures.  Note,  this  testing 
accepted  a  wide  margin  for  error  due  to  the  relatively 
unsophisticated  fracture  mechanics  methods  available  for 
analyzing  components  subjected  to  high  temperature  and  low 
frequency  cycle  loading  [3] . 

The  life  cycle  of  a  turbine  disk  is  a  complex  one 
consisting  of  frequent  load  cycle  variations  with  ambient 
temperatures  of  up  to  1350°  F.  The  elevated  temperatures 
introduce  time  dependent  creep  phenomenon  which  interact 
with  varied  load  spectra  to  produce  complex  material 
behavior  [4] .  If  the  material  behavior  could  be  determined 
for  the  typical  life  cycle  then  accurate  remaining  life 


CRACK  CLOSURE  CHARACTERISTICS  CONSIDERING 
CENTER  CRACKED  AND  COMPACT 
TENSION  SPECIMENS 


I .  Introduction 

Background 

The  growing  use  of  expensive  high  performance  gas 
turbine  engines  in  multi-million  dollar  aircraft  has 
created  a  problem  in  the  United  States  Air  Force.  An 
accurate  failure  prediction  method  is  unavailable  to  retire 
critical  high  temperature  jet  engine  components  [1] . 

Normally,  aircraft  components  are  periodically 
inspected  for  flaws  and  returned  to  service  if  the  flaws 
can  not  grow  to  critical  size  prior  to  subsequent  periodic 
inspections.  However,  critical  engine  components,  like 
turbine  engine  disks,  are  removed  from  service  at  a  time 
when  statistically  1  in  1,000  would  be  expected  to  initiate 
a  flaw  of  some  finite  length  (0.03  in.).  From  a  safety 
standpoint,  this  policy  works  well.  However,  by  definition 
99.9  percent  of  the  retired  disks  still  have  useful  life. 
Over  80  percent  of  the  disks  have  at  least  ten  lifetimes 
remaining  and  over  50  percent  have  at  least  25  lifetimes 
remaining  [1]  . 

If  one  could  accurately  predict  the  remaining  life  of  a 
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Abstract 


Due  to  growing  use  of  expensive,  high  performance  gas 
turbine  engines  in  the  United  States  Air  Force,  there  is  a 
need  for  improved  failure  prediction  methods  for  critical 
high  temperature  engine  components.  This  new  study  expands 
current  research  in  the  area  of  high  temperature,  low  cycle 
fatigue  of  IN-100  at  1350°  F,  the  superplastic  alloy  used 
in  F-100  engine  turbine  disks. 

An  in-house,  2-D,  finite  element  program  named  VISCO 
employs  the  Bodner-Partom  Constitutive  equation  to 
accurately  model  the  principal  features  of  completely 
reversed  cyclic  loading.  VISCO  is  used  to  compare  the 
effects  on  material  behavior  by  considering  a  2.5  Hz. 
compact  tension  specimen,  a  .167  Hz.  compact  tension 
specimen,  and  a  2.5  Hz.  center  cracked  specimen  subjected 
to  fully  reversed  cyclic  loading  with  a  stress  intensity 
factor  of  35  and  45  ksi  /  in.  The  comparisons  point  out  the 
findings  of  Linear  Elastic  Fracture  Mechanics  must  be 
modified  under  conditions  of  high  temperature 
viscoplasticity. 


£  o i i  during  a  given  time  step.  Stress/ strain  tolerances 
effect  time  stepping  as  follows: 


P_  = 


a  1  -a1"1 
e _ e 

„  i-1  _ 

ae  TOL 


P  = 


(d£  P)1 


eTOTAL  £T0L 


(3.1) 


(3.2) 


where  superscript  i  refers  to  the  timestep,  is 
effective  stress,  is  effective  plastic  strain  and 

£TOL  defined  as 


£  TOTAL 


+  o.5  e  2y 

xy 


(3.3) 


PQ  and  P£  are  parameters  evaluated  for  each  element  and  the 
largest  of  the  two  is  set  equal  to  P.  The  timestep  size  is 
then  determined  from  the  following  equations  [8]: 


dt1  =  0.8  dt1/ P 

if 

P  >  1 

dt1  =  dt1'1 

if 

0.8  <  P  <  1 

dt1  =  1.25  dt1"1 

if 

0.65  <.  P  < 

dt1  =  1.5  dt1"1 

if 

P  <  0.65 

thus,  the  timestep  interval  is  increased  or  decreased  for 
efficient  and  accurate  computations. 


To  improve  understanding  of  the  material  behavior  of 
critical  high  temperature  engine  components,  the  full  range 


of  possible  fatigue  load  spectra  must  be  represented. 
Consequently,  a  completely  reversed  (R-ratio  =  -1.0) 
sawtoothed  load-time  pattern  is  used  to  drive  the  finite 
element  models.  Frequency,  max  load,  and  R-ratio  can  be 
individually  prescribed  for  each  computer  simulation.  A 
typical  load  spectrum  is  shown  in  Figure  3.1. 


„  . .  min  load  ,  „ 

R-ratio  =  - : — v  =  -1.0 

max  load 

Figure  3.1  Typical  Load  Cycles  2.5  Hz, 


Other  than  minor  modifications  listed  in  Appendix  A, 
VISCO,  as  used  by  Wilson  [11]  was  unchanged. 


Finite  Element  Modeling 

Two  finite  element  models  shown  in  Figures  3.2  and  3.3 
were  used  in  this  study:  the  382  element  compact  tension 
specimen  (CTS)  and  the  355  element  center  cracked  specimen 
used  by  Hinnerichs  [6] .  Due  to  symmetry  only  the  top  half 
of  the  CTS  was  modeled.  Adjacent  elements  differ  in  size 
by  no  more  than  a  factor  of  2. 


eg 


aeasaai 


gUjpg 

mmM 


Figure  3.2  382  Element  Compact  Tension  Specimen  Mesh 
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2.5  IN 
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Figure  3.3  355  Element  Center  Cracked  Specimen  Mesh 


The  fine  mesh  elements  near  the  crack  tip  have  an  area  of 
—  6  2 

4.8848  x  10  in.  .  The  crack  length  and  specimen 
thickness  were  .6630  in.  and  .2154  in.  respectively.  Only 
the  top  right  quarter  of  the  center-cracked  specimen  was 
modeled  due  to  symmetry.  Like  the  CTS,  adjacent  elements 
differed  in  size  by  no  more  than  a  factor  of  2.  The  fine 
mesh  elements  near  the  crack  tip  were  much  smaller  than  the 
CTS  with  an  area  of  3.051757  x  10  ^  in.^.  The  crack 
length,  2a,  is  .2734  in.  and  the  plate  thickness  is  .3000 
in . 


The  loads  were  adjusted  on  each  model  to  provide  the 
stress  intensity  factor  ( K ^ )  as  shown  in  Table  3.1.  Kj. 
for  the  CTS  was  found  using 


2  a  3 

t  29.6  (f)  -  185.5  (^)  +  655.7  (J) 


4  5 

-  1017  (^)  +  638.9  (f)  ] 


(3.4) 


where  P  is  load,  B  is  thickness,  a  is  crack  length  and  w  is 
width,  [18].  Kj  for  the  center  cracked  specimen  was 
found  using 

K,  =  y  y~na  (sec  — )  (3.5) 

I  A  w 


where  A  is  the  area  over  which  the  load  is  applied  [18]. 
See  Figures  3.4  and  3.5  for  dimensions. 
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Figure  3.4  Compact  Tension  Specimen  Geometry 


K (ksi/ 1 In . ) 


CTS  Load  (Lbs.) 


Center  Cracked 
Plate  Load  (lbs.) 


35  1710.0 

45  2198.0 


8029.9 

10300.0 


Table  3.1  Specimen  Loads  and  Stress  Intensity  Factors 


A  stress/strain  tolerance  investigation  was  performed 
to  find  a  balance  between  accuracy  and  computer  computation 
time.  Several  analyses  under  plane  stress  conditions  were 
performed  using  different  values  of  stress  and  strain 
tolerances.  The  results  from  both  meshes  repeated 
Hinnerichs'  findings  that  computer  time  increases  rapidly 
as  the  stress/strain  tolerances  are  reduced  to  .01  [6] .  It 
was  also  noted  that  the  change  in  plastic  work,  from  one 
tolerance  analysis  to  the  next,  for  each  of  the  cases  in¬ 
vestigated,  diminished  as  .01  was  approached,  See  Figure 
3.6.  One  may  observe  from  Figure  3.6  that  the  center 
cracked  specimen  is  least  affected  by  the  change  in  stress/ 
strain  tolerances.  This  can  be  attributed  to  smaller  area 
per  element  near  the  crack  tip.  Using  stress  and  strain 
tolerances  of  .01  as  a  datum  point,  tolerances  were  in¬ 
creased  to  provide  no  more  than  a  2%  error  in  plastic  work 
over  one  computer  loading  cycle,  this  allowed  a  minimum  of 
5  load  cycles  in  2500  seconds  of  central  processor  time  on 
the  CDC  CYBER  845  computer  for  all  cases  studied.  This 
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PERCENT  ERROR 


FIG  3.6  ERROR  VS.  STRESS/STRAIN  TOLERANCE  SETTING 
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process  established  the  tolerances  for  the  CTS  and  the 
center  cracked  specimen.  Stress/ strain  tolerances  in  this 
study  were  set  at  .06  and  .10  for  the  CTS  and  center 
cracked  specimen  respectively. 
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IV.  Results  and  Discussion 


Throughout  this  chapter,  comparisons  will  be  made  of 
three  cases  under  two  different  stress  intensity  factors, 
35  and  45  ksi  /  in : 

(1)  CTS  with  2.5  Hz.  Load  Cycling 

(2)  CTS  with  .167  Hz.  Load  Cycling 

(3)  Center  cracked  specimen  with  2.5  Hz.  Load 
Cycling 

All  three  cases  model  the  material  properties  of  IN-100  at 
1350°  F  subjected  to  a  completely  reversed  (R-ratio  = 

-1.0)  cyclic  load.  Previous  work  by  Wilson  [11],  which 
dealt  with  only  case  (1),  was  verified  as  part  of  this 
study.  The  desire  to  investigate  interesting  findings  in 
Wilson's  work  with  the  2.5  Hz.  CTS,  along  with  the  con¬ 
tinuing  need  for  better  understanding  of  high  temperature, 
low  cycle  fatigue,  led  to  the  inclusion  of  the  other  two 
cases.  As  discussed  in  Chapter  III,  steps  were  taken  to 
insure  the  same  degree  of  accuracy  and  same  stress  inten¬ 
sity  factors  exists  in  each  case  so  valid  comparisons  can 
be  made.  The  results  are  separated  into  three  areas: 

(a)  Crack  opening  displacements  behind  the  crack 
tip  at  full  negative  load. 

(b)  Profiles  of  y-stress  and  y-strain  fields  in 
front  of  the  crack  tip,  (the  y  components 
along  the  horizontal  line  of  symmetry  were 


chosen  since  they  are  the  most  significant) , 


and 

(c)  Plastic  zone  size  and  shape  estimations. 

Closure  Behind  Crack  Tip  at  Full  Negative  Load 

One  interesting  finding  in  Wilson's  CTS  work  [11] 
considered  the  fact  that  an  incomplete  closure  behind  the 
crack  tip  (that  is,  the  open  side)  existed  at  full  negative 
load.  Looking  at  the  way  the  CTS  is  loaded  (Figure  3.3), 
it  can  be  seen  that  no  load  is  applied  directly  over  the 
near  field  of  the  crack  tip.  In  fact,  the  loading  geometry 
is  much  like  a  cantilever  beam.  Wilson  proposed  that  the 
plastic  deformation  near  the  crack  tip  from  the  full 
positive  load,  acts  like  a  fulcrum  and  prevents  full 
closure.  Further,  it  seems  reasonable  to  assume  that  a 
center  cracked  specimen,  with  uniform  loading  directly  over 
the  near  field  of  the  crack  tip  (Figure  3.4),  would  close 
completely  upon  full  negative  loading.  Figures  4.1  through 
4.4  refute  this  idea.  Even  though  the  2.5  Hz.  center 
cracked  specimen  is  characterized  by  more  full  closure  than 
the  CTS  at  2.5  Hz.  and  .167  Hz.,  the  highly 
viscoplastical ly  deformed  area  surrounding  the  crack  tip 
still  prevents  full  closure.  The  .167  Hz.  CTS  has  the 
least  closure  and  the  2.5  Hz.  CTS  closure  is  in  between. 

One  should  notice  the  difference  that  a  number  of  cycles 
makes.  This  change  can  be  observed  in  Figures  4.1  and  4.2 


I 


i 

i 


(or  4.3  and  4.4  at  45  ksi/  in);  the  lack  of  full  closure 
on  the  first  full  negative  load  (.75  cycles)  becomes 
slightly  worse  as  cycles  increase  to  the  fourth  full 
negative  load  (3.75  cycles).  Also  note  the  difference 
caused  by  the  stress  intensity  factors  in  Figure  4.1  versus 
Figure  4.3  and  Figure  4.2  versus  Figure  4.4.  The  lack  of 
full  closure  increases  significantly  as  stress  intensity 
factor  is  increased  from  35  to  45  ksi/  in.  The  reasons  for 
this  phenomenon  becomes  evident  by  examining  the  state  of 
stress  and  strain  in  front  of  the  crack  tip. 

y-Stress  and  Total  y-Strain  Fields  in  Front  of  Crack  Tip 

The  y-stress  versus  total  y-strain  curves  for  all  three 
cases  at  both  35  and  45  ksi/  in  are  presented  in  Figures 
4.5  through  4.10.  These  figures  plot  the  stress/strain 
values  in  the  near  field  of  the  crack  tip  (.004  inches  in 
front  of  the  crack  tip).  Comparing  Figure  4.5  with  Figure 
4.6,  one  can  see  the  effect  of  varying  only  cyclic 
frequency.  At  2.5  Hz.,  the  full  load  has  0.1  seconds  to  be 
applied  while  the  .167  Hz.  case  has  1.5  seconds.  The 
result  is  higher  positive  stress  and  lower  total  strain  in 
the  2.5  Hz.  case.  This  phenomenon  is  known  as  rate  sensi¬ 
tivity.  The  .167  Hz.  specimen,  consequently,  experiences 
more  plastic  strain  which  is  the  reason  why  there  is  less 
closure  behind  the  crack  tip  at  full  negative  load. 

This  larger  region  of  plastic  strain  results  in  greater 
compressive  stresses  (Figure  4.6)  as  the  surrounding 
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material  attempts  return  to  its  pre-strained  state. 

Looking  at  Figure  4.7,  the  center  cracked  specimen  develops 
approximately  the  same  maximum  tensile  stress  as  the  2.5 
Hz.  CTS ,  but  the  total  y-strain  and  the  maximum  y-com- 
pressive  stress  is  much  greater  than  either  CTS  cases.  The 
uniform  distributed  loading  in  the  positive  direction 
acting  on  the  center  cracked  specimen  increases  the  tensile 
y-stress  field  throughout  the  material  and  causes  much  more 
plastic  strain  than  either  CTS  case.  At  this  point,  one 
would  expect  less  closure  behind  the  crack  tip,  at  full 
negative  load,  than  even  the  .167  Hz.  CTS.  However,  the 
uniform  distributed  load  in  the  negative  direction  in¬ 
creases  the  compressive  y-stress  field  throughout  the 
material.  This  compressive  y-stress  field  does  more  work 
in  overcoming  the  larger  plastic  deformation  to  produce  the 
most  closure  behind  the  crack  tip.  The  large  compressive 
stresses  in  the  2.5  Hz.  center  cracked  specimen,  depicted 
in  Figure  4.7,  are  further  evidence  of  this  phenomenum.  Of 
the  three  cases,  the  center  cracked  specimen  experiences 
the  widest  range  of  stresses  and  total  strains  in  the  near 
field  of  the  crack  tip. 

Increasing  the  stress  intensity  factor  from  35  to  45 
ksi/Tn,  as  shown  in  Figures  4.8  -  4.10,  has  the  expected 
effect  of  increasing  tensile  and  compressive  stresses; 
which  in  addition,  produces  a  greater  amount  of  comp¬ 
ressive  plastic  strain.  It  is  this  greater  amount  of 
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FIG  4.1  Y-DI5PL  BEHINO  CRACK  (FULL  NEG.  LOAD) 
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compared  with  the  region  of  negative  plastic  strain  under 
full  negative  load,  at  35  ksi/ in.  Figures  4.33  through 
4.35  represent  the  45  ksi/ in  cases.  The  region  of 
negative  plastic  straining  is  approximately  one  fourth  the 
area  of  the  region  of  positive  straining  in  all  cases. 

These  regions  can  be  thought  of  as  residual  compressive 
stress  zones,  since  the  residual  stress  must  be  overcome 
before  the  region  can  be  plastically  strained  in  tension. 
The  result  of  these  residual  compressive  stress  zones,  when 
large  enough,  is  the  redistribution  of  stresses  in  the  form 
of  the  discontinuities  seen  in  the  y-stress  versus  distance 
in  front  of  the  crack  tip  figures.  Note  that  the  region  of 
negative  plastic  straining  in  the  two  2.5  Hz.  center 
cracked  specimen  cases  and  the  45  ksi  /  in  CTS  at  .167  Hz. 
case  (Figures  4.32,  4.35,  and  4.34)  have  the  size  and 
position  corresponding  to  the  size  and  position  of  the 
y-stress  discontinuities  in  Figures  4.13  and  4.14. 
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the  raw  data  printout  of  Bodner-Partom  plastic  strain.  It 
is  noteworthy  that  the  region  closest  to  the  crack  tip  does 
not  return  to  a  state  of  zero  strain  even  during  full 
negative  load. 

One  way  to  define  compressive  plastic  straining  is  to 
compare  total  strain  at  full  negative  load  with  total 
strain  at  full  positive  load.  When  the  difference  exceeds 
the  appropriate  value  in  Table  4.2,  compressive  plastic 
straining  has  occurred.  Unlike  Table  4.1,  Table  4.2 
accounts  for  the  full  range  of  elastic  strain,  that  is, 
elastic  strain  in  tension  plus  elastic  strain  in  com¬ 
pression.  The  results  in  Table  4.2  can  be  verified  by 
estimating  the  linear  range  of  the  stress  strain  curves 
Figures  4.5  through  4.10. 


Case 

Plastic 

Strain 

in/in 

Total  Strain 
in/in 

35  ksi  /  in 

Total  Strain 
in/in 

45  ksi  /  in 

(1) 

2.5  Hz.  CTS 

.001 

.013 

.013 

(2) 

.167  Hz.  CTS 

.001 

.011 

.011 

(3) 

2.5  Hz.  Ctr. 

Crk.  .001 

.013 

.013 

Table  4.2  Total  Elastic  Strain  Range  from  Full  Tension 

to  Full  Compression 

Plastic  Zone 

The  actual  size  and  shape  of  the  plastic  zone  will  help 
explain  the  material  behavior.  Figures  4.30  through  4.32 
show  the  regions  of  plastic  strain  under  full  positive  load 


Bodner-Partom  results,  will  be  used  to  define  the  yield 
point  and  the  beginning  of  the  plastic  zone.  Table  4.1 
lists  the  total  strain  occurring  in  the  three  cases  when 
0.1%  plastic  strain  has  been  achieved.  These  results  can 
be  verified  by  reviewing  the  y-stress/y-strain  curves 
(Figures  4.5  through  4.10). 


Case 

Plastic 

Strain 

in/in 

Total 

Strain 

35  ksi/TrT 

Total 

Strain 

45  ksi/Th 

(1) 

2.5  Hz.  CTS 

.001 

.007 

.007 

(2) 

.167  Hz.  CTS 

.001 

.006 

.006 

(3) 

2.5  Hz.  Ctr.  Crk. 

.001 

.007 

.007 

Table  4.1  Total  Elastic  Strain  Range  Under  Tension 

The  .167  Hz.  Case,  with  the  .001  plastic  strain 
removed,  matches  the  normal  elastic  range  of  .5%  based  on 
observations  in  the  theory  of  elasticity  [5] .  The  higher 
values  at  2.5  Hz.  are  due  to  the  rate  sensitivity  of  the 
material;  that  is,  high  loading  rates  push  up  the  yield 
point  since  less  time  is  available  for  plastic  straining  to 
redistribute  the  stress. 

Using  the  values  in  Table  4.1,  an  approximation  can  be 
made  to  determine  where  each  strain  field  goes  plastic 
under  tension  by  using  the  full  positive  load  portion  of 
Figures  4.28  and  4.29.  This  method  is  a  good  approximation 
for  full  positive  load  plastic  straining  when  compared  with 
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discontinuity.  Further  investigation  considering  strain  in 
front  of  the  crack  is  warranted. 

The  total  y-strain  field  at  4.25  cycles  (full  positive 
load)  for  all  three  cases  is  compared  with  the  total 
y-strain  field  at  3.75  cycles  (full  negative  load),  in 
Figure  4.28  for  35  ksi  /  in  and  Figure  4.29  for  45  ksi  /  in. 
The  total  y-strain  appears  to  increase  exponentially  close 
to  the  crack  tip.  This  corresponds  to  the  exponential  rise 
in  stresses  close  to  the  crack  tip  predicted  by  Linear 
Elastic  Fracture  Mechanics.  The  Bodner-Partom  Constitutive 
Law,  with  its  ability  to  model  viscoplastic  material 
behavior,  allows  the  prediction  of  this  phenomenon.  The 
center  cracked  specimenexhibits  the  highest  y-strain 
throughout  the  near  field  of  the  crack  tip  under  full 
positive  load,  and,  except  for  the  region  right  next  to  the 
crack  tip  in  the  35  ksi/  in  case,  the  center  cracked 
specimen  has  the  lowest  value  of  strain  at  full  negative 
load.  The  load  geometry  for  this  case  puts  the  center 
cracked  specimen  through  a  much  wider  range  of  strains  than 
the  CTS  loading  geometry. 

A  criteria  is  needed  to  determine  when  plastic 
straining  occurs.  Conventional  strength  of  materials 
normally  defines  yield  point  as  .2%  permanent  deformation 
represented  in  a  uniaxial  stress/ strain  curve  [19] .  Con¬ 
sequently,  the  plastic  strain  components  in  the  y  direction 
will  be  used  to  predict  when  an  element  undergoes  signif¬ 
icant  plastic  strain.  0.1%  plastic  strain,  based  on  the 


discontinuity.  To  better  illustrate  the  effect  of  cycling 
on  the  y-stress  field  near  the  crack,  tip,  the  y-stress  at 
the  first  positive  load  is  compared  with  the  y-stress  for 
the  fifth  full  positive  load  for  each  case  at  35  and  45 
ksi/irT.  -  see  Figures  4.20  through  4.25.  These  six  figures 
clearly  demonstrate  the  viscoplastic  property  of  stress 
redistribution  due  to  time  dependent  plastic  straining. 

The  highest  stresses  near  the  crack  tip  are  gradually 
redistributed  through  time-dependent  plastic  strain.  For 
example,  in  Figure  4.10  the  stress  field  at  4.25  cycles 
starts  out  with  lower  stress  values,  but  drops  off  less 
rapidly  than  the  .25  cycle  results. 

The  stress  field  from  the  fourth  full  negative  load 
(3.75  cycles)  is  superimposed  on  the  stress  field  from  the 
fifth  full  positive  load  (4.25  cycles)  in  Figures  4.26  and 
4.27.  This  way,  it  will  be  easier  to  see  the  effect  of  a 
negative  load  stress  field  on  the  subsequent  positive  load 
stress  field.  As  expected,  the  discontinuities  coincide 
with  the  maximum  compression  stresses.  The  center  cracked 
specimen  at  35  ksi  / in  and  45  ksi  /in  along  with  the  .167 
CTS  at  45  ksi  / in  have  higher  magnitudes  of  compressive 
stress  than  the  other  cases  without  discontinuities. 
Presumably,  these  stresses  are  large  enough  to  cause 
negative  plastic  deformation  with  resulting  residual 
compressive  stresses.  These  residual  stresses  effectively 
subtract  from  the  maximum  tensile  stresses  forming  the 


Kj  does  not  produce  the  same  stress  near  the  crack  tip 
when  applied  at  different  cyclic  frequencies;  now  it  is 
evident  that  the  same  Kj  does  not  produce  the  same  stress 
near  the  crack  tip  when  applied  to  different  load 
geometries.  Now  that  the  stress  field  in  the  so-called 
virgin  material  has  been  analyzed,  what  is  the  effect  of 
cycling  with  completely  reversed  load  on  the  stress  field? 

After  approximately  five  complete  cycles  at  35  ksi/ in, 
the  y-stress  field  under  full  positive  load  has  gone 
through  some  changes.  Looking  at  Figure  4.13,  one 
immediately  notices  a  discontinuity,  or  sudden  dip,  in  the 
center  cracked  specimen.  The  discontinuity,  .004  inches 
from  the  crack  tip,  was  not  expected.  Figure  4.14,  at  45 
ksi/Tn,  shows  the  discontinuity  has  become  more  pronounced 
and  has  moved  further  away  from  the  crack  tip  (.008 
inches).  The. 167  Hz.  CTS  seems  to  be  developing  a  discon¬ 
tinuity  also.  Note  also.  Figures  4.13  and  4.14  verify  the 
conclusion  that  the  same  does  not  produce  the  same 
y-stress  near  the  crack  tip  when  applied  to  different  load 
geometries.  Figures  4.15  through  4.19  trace  the  development 
of  the  center  cracked  specimen  discontinuity  over  5  cycles 
at  45  ksi  /  in.  Clearly  the  size  of  the  discontinuity  is 
related  to  the  magnitude  of  the  compressive  stress  at  full 
negative  load.  It  will  be  shown  that  negative  .  stic 
straining,  produced  by  the  compressive  stress,  results  in  a 
region  of  residual  compressive  stress  corresponding  to  the 


considering  R-ratio  =  -  1.0,  results  close  to  the  fifth 
cycle  will  be  considered  stabilized.  The  next  series  of 
figures  allow  a  closer  examination  of  the  stress  and  strain 
field  in  front  of  the  crack  tip. 

Figure  4.11  shows  y-stress  as  a  function  of  distance  in 
front  of  the  crack  tip  at  full  positive  load  on  the  first 
cycle.  The  center  cracked  specimen  curve  verifies  that  a 
higher  stress  field  exists  throughout  the  specimen  compared 
to  both  CTS  curves.  Note  however,  the  2.5  Hz.  CTS 
approaches  convergence  with  the  center  cracked  specimen 
curve  close  to  the  crack  tip.  This  seems  to  verify  that 
the  same  Kj  will  produce  similar  characteristics  in 
different  structures.  The  contrary  of  this  will  become 
apparent  in  the  subsequent  analysis  of  this  study's 
results.  The  2.5  Hz.  CTS  has  the  expected  higher  stress 
level  compared  to  the  .167  Hz.  CTS  close  to  the  crack 
tip.  This  results  from  the  rate  sensitivity  discussed 
earlier.  The  2.5  Hz.  CTS  curve  drops  rapidly  below  the 
.167  Hz.  CTS  curve  and  rejoins  away  from  the  crack  tip. 

The  explanation  for  this  lies  behind  the  lower  load  rate 
which  allows  more  time  for  the  plastic  strain  process  in 
the  .167  Hz.  CTS  to  distribute  the  stress  over  more 
elements,  this  process  produces  lower  stress  close  to  the 
crack  tip  and  slightly  higher  stress  in  neighboring 
elements.  Figure  4.12,  at  45  ksi/TrT,  essentially  magnifies 
the  35  ksi/ in  results.  It  is  already  clear  that  the  same 


compressive  plastic  strain  that  results  in  even  less 
closure  behind  the  crack  tip  than  in  the  35  ksi/in  cases. 
Notice  now,  especially  in  Figures  4.9  and  4.10,  that 
negative  plastic  straining  is  occurring  in  the  compressive 
region  of  the  stress/strain  curves.  This  is  causing  an 
overlapping  of  the  curves  as  each  cycle  enters  the  region 
of  full  negative  load.  The  ability  to  account  for  plastic 
straining  in  compression  as  well  as  plastic  strain  in 
tension  is  a  characteristic  of  the  Bodner  material 
idealization  [10] . 

Reviewing,  in  order.  Figures  4.5  through  4.10,  the 
maximum  compressive  stresses  progressively  approach  the 
maximum  tensile  stresses  in  magnitude,  this  is  a  direct 
result  of  the  increasing  size  of  the  plastic  zone,  resists 
closure  and  the  return  to  the  undeformed  state.  Finally, 
compressive  stresses  have  grown  large  enough,  especially  in 
Figures  4.9  and  4.10,  to  cause  plastic  strain  in  the 
negative  direction. 

In  all  six  stress/strain  curves  there  is  a  noticeable 
ratcheting  (stepped  movement)  to  the  right.  The  rate  of 
change  in  the  ratcheting  diminishes  as  cycles  increase  to  a 
point  where  further  plastic  straining  may  cease.  Work  by 
Wilson  [11]  showed  that,  for  an  R-ratio  of  0.1,  this  stable 
zone  occurred  after  23  cycles.  Due  to  the  prohibitive 
central  processor  time,  associated  with  completing  the 
possible  high  number  of  cycles  for  a  stable  plastic  zone 
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FIG  4.11  Y-STRESS  IN  FRONT  OF  CROCK  (FULL  POS  LD  ) 
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FIG  4.15  Y-STRESS  IN  FRONT  OF  CRACK 
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FIG  4.19  Y-STRESS  IN  FRONT  OF  CRRCK 
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FIG  4.32  REGION  OF  PLASTIC  STRAINING 


FIG  4.30  REGION  OF  PLASTIC  STRAINING 


65 


FIG.  4.27  Y-STRESS  IN  FRONT  OF  CRACK  TIP 
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FIG  4.35  REGION  OF  PLRSTIC  STROKING 
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V.  Conclusions 


The  following  conclusions  for  IN-100  at  1350°  F  can 
be  made  relative  to  the  work  carried  on  within  the  study: 

1.  None  of  the  three  cases,  at  either  35  or  45 
k si/in,  close  completely  in  the  region  immediately  behind 
the  crack  tip,  even  at  full  negative  load.  This  phenomenon 
is  a  result  of  the  large  region  of  plastic  strain  formed 
during  full  positive  load. 

2.  As  the  number  of  cycles  increase  and  the 
region  of  plastic  strain  fully  develops,  slightly  less 
closure  occurs  behind  the  crack  tip. 

3.  As  Kj  is  increased,  the  region  of  plastic 
strain  grows,  and  signif icantly  less  closure  occurs 
behind  the  crack  tip. 

4.  The  amount  of  closure  for  the  center  cracked 
specimen  is  greater  than  both  CTS  cases,  but  significantly, 
it  does  not  close  all  the  way. 

5.  The  size  of  the  plastically  strained  region  at 
the  crack  tip  is  a  major  factor  determining  the  amount  of 
closure  behind  the  crack  tip.  Load  and  specimen  geometry, 
in  the  particular  case  of  the  center  cracked  specimen  as 
compared  with  both  CTS  cases,  has  an  overriding  influence 
on  closure.  Furthermore,  due  to  the  uniform  loading  in  the 
center  cracked  specimen,  it  experienced  a  wider  range  of 
stresses  and  strains  in  both  compression  and  tension 
throughout  the  material. 
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6.  Maximum  stresses  in  front  of  the  crack  tip 
decrease  and  approximately  stabilize  at  a  lower  value  after 
5  cycles.  This  is  a  result  of  the  time  dependent  plastic 
strain  and  the  redistribution  of  stresses. 

7.  The  same  stress  intensity  factor,  Kj,  does 
not  produce  the  same  stresses  near  the  crack  tip  when 
applied  to  different  geometries  of  frequencies.  The 
differences  increase  as  Kj  increases  and/or  as  the  number 
of  cycles  increases. 

8.  Lower  maximum  stresses  and  increased  total 
strain  in  front  of  the  crack  tip  are  achieved  when  a  lower 
cyclic  load  frequency  is  used.  This  is  the  characteristic 
of  load  rate  sensitivity,  which  is  due  to  the  redistribu¬ 
tion  of  stresses  through  time-dependent  plastic  strain. 

9.  After  cycling  at  a  -1.0  R-ratio,  a  discon¬ 
tinuity,  or  small  region  of  significantly  lower  stress  in 
the  stress  field  in  front  of  the  crack  tip,  develops.  This 
is  a  result  of  residual  compressive  stresses  at  the  discon¬ 
tinuity  location  caused  by  negative  plastic  straining 
during  full  negative  load.  This  phenomenon  occurred  only 
in  the  center  cracked  specimen  at  35  and  45  ksi/in  and  the 
45  ksi/in"  CTS  at  .167  Hz.  Also,  these  cases  exhibited  the 
most  amount  of  negative  plastic  straining  during  full 
negative  load. 

10.  An  analysis  of  the  region  undergoing  plastic 
straining  near  the  crack  tip,  at  full  positive  and  full 


negative  loads,  verifies  that  the  discontinuity  is  a  result 
of  residual  compressive  stress. 

11.  The  region  of  compressive  plastic  straining 
during  full  negative  load  is  approximately  one  fourth  the 
area  of  the  region  tensile  plastic  straining  during  full 
positive  load. 
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Appendix  A 


The  computer  program,  VISCO,  modified  by  Wilson  [11] 
for  negative  R-ratio  cyclic  loading  was  changed  in  several 
areas.  The  changes,  while  significant  to  the  results  of 
this  study,  did  not  interfere  with  the  proven  Bodner-Partom 
viscoplastic  computation  procedure.  This  fact  was  verified 
by  duplicating  Wilson's  results  using  this  version  of 
VISCO. 

The  modifications  to  VISCO  in  this  study  can  be  divided 
into  three  major  categories:  cosmetics,  specimen 
adaptation  and  output  enhancement.  Cosmetics  included 
commenting  out  unused  code  (00276-00287,  00299,  00300, 
00908-00913)  and  removal  of  unused  subroutines  (geom,  func 
and  disp) .  The  category  of  specimen  adaptation  deals  with 
the  integral  boundary  conditions  originally  inserted  by 
Wilson  [11]  to  prevent  negative  displacements  along  the 
crack  edge  during  negative  cycling  (00697-00700) .  This 
modification  was  necessary  to  represent  the  center  cracked 
specimen  boundary  conditions.  Further,  the  center  cracked 
specimen  mesh  had  different  elements  of  interest  in  front 
of  the  crack  tip  (00864-00865,  00976-00994,  01071-01072). 
Output  enhancement  included  format  changes  (01003,  01006, 
01008)  and  output  limitation  (00936-00939,  00973-00974) 
which  combined  to  make  data  more  accessible  for 
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analysis.  One  modification  added  Z,  the  Bodner-Partom 
hardness  parameter  to  the  output  (00866) .  The  final 
modification  output  enhancement  turned  out  to  be  the  most 
significant . 

Due  to  the  nature  of  this  study,  result  comparison  of 
unique  cases,  it  was  critical  that  output  of  results  could 
be  directed  consistently  at  the  precise  moment  of 
interest.  The  original  timestep  algorithm  used  to  direct 
output  had  a  built-in  inaccuracy  which  accrued  each  time 
results  were  printed.  For  example,  when  the  moment  of 
interest  was  full  positive  load  after  two  cycles,  data  at 
+83.7%  load  was  printed.  Further,  the  subsequent  full 
negative  load  data  was  actually  based  on  only  -71.5%  load. 
Since  this  error  accrued  as  cycles  increased,  output  close 
to  desired  cyclic  load  was  hard  to  come  by.  After  two 
simple  modifications:  adding  line  00253  and  changing  line 
00291  from  "TP=0",  to  the  present  form,  loads,  even  in  the 
worst  case  (.167  Hz.  CTS) ,  were  never  off  more  than  4%  of 


desired  load. 
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